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The magnetic state of a quantum dot attached to superconducting leads is experimentally shown 
to be controlled by the superconducting phase difference across the dot. This is done by probing 
the relation between the Josephson current and the superconducting phase difference of a carbon 
nanotube junction whose Kondo energy and superconducting gap are of comparable size. It exhibits 
distinctively anharmonic behavior, revealing a phase mediated singlet to doublet transition. We 
obtain an excellent quantitative agreement with numerically exact quantum Monte Carlo calcula¬ 
tions. This provides strong support that we indeed observed the finite temperature signatures of the 
phase controlled zero temperature level-crossing transition originating from strong local electronic 
correlations. 
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When a localized magnetic moment interacts with a 
Fermi sea of conduction electrons, the Kondo effect can 
develop: spin-flip processes lead to a many-body singlet 
state in which the delocalized electrons screen the mo¬ 
ment. Quantum dots (QD) in the Coulomb blockade 
regime and particularly carbon nanotube (CNT) dots 
constitute ideal systems for the investigation of Kondo 
physics at the single spin level nHa. In these systems, it 
is possible to control the number of electrons on the dot 
varying a gate voltage. For an odd occupation, the dot 
accommodates a magnetic moment which is screened pro¬ 
vided that this is not prohibited by an energy scale larger 
than the Kondo energy Temperature is the most 

obvious obstacle to the development of the Kondo effect 
since Tx can be smaller than IK. However, if temper¬ 
ature is sufficiently low, the Kondo effect may compete 
with other quantum many-body phenomena such as su¬ 
perconductivity, for which the formation of Cooper pairs 
of energy A may prevent the screening of the dot’s spin. 
This situation can be investigated using superconducting 
hybrid junctions, where a supercurrent is induced by the 
proximity effect, for example in CNT-based QDs |4] or 
semiconductor-based ones [5). 

A setup of a high resistance tunnel barrier between 
two superconductors, also called Josephson junction (JJ), 
carries a supercurrent I = Ic sin cp, with the critical cur¬ 
rent Ic- The superconducting phase difference across 
the junction ip controls the amplitude and the sign of the 
supercurrent. This is the Josephson relation, the most 
famous example of a current-phase relation (CPR). In 
some peculiar systems such as ferromagnetic supercon¬ 
ducting junctions, the transmission of Cooper pairs gives 
rise to a TT phase shift of the CPR [6]. In QD JJs (tunnel 
barrier replaced by QD) in the strong Coulomb blockade 
regime where the Kondo effect is negligible, such a tt shift 
is observed as well since the tunneling of a Cooper pair 


implies reversing the order of particles within this pair. 
This leads to a gate-controlled sign reversal of the CPR 
when the parity of the number of electrons is changed, as 
was observed experimentally IMHl- In contrast, if the 
Kondo effect and thus local correlations prevail, the spin 
of the dot is screened by unpaired electrons leading to 
a singlet ground state: the 0-junction is then recovered 
even though the parity of the dot is still odd. 

The switching from 0- to 7r-junction behavior as a 
function of a variety of energy scales of the QD JJ in 
the presence of local correlations was extensively stud¬ 
ied theoretically [HHini. The scales are the broadening 
F of the energy levels in the dot due to the coupling to 
the reservoirs, the superconducting gap A of the con¬ 
tacts, the dot’s charging energy U and its level energy 
e. When these parameters fall into the O-tt transition 
regime, it was predicted that the ground state of the 
system—singlet or doublet—depends on the phase dif¬ 
ference (p, undergoing a level-crossing transition. This 
leads to a characteristic anharmonicity of the CPR for 
temperature T > 0 and a jump at a critical phase pc for 
T = 0. In other words, in this particular regime of pa¬ 
rameters, the magnetic state of the dot is predicted to be 
governed by the superconducting phase difference across 
the junction. 

This Kondo related O-tt transition was earlier observed 
experimentally as a function of the gate voltage [3 El [22] 
and the spectroscopy of Andreev bound states enabled 
a better understanding of the involved physics [23H27] . 
Measurements of the CPR of a QD JJ embedded in a 
SQUID were also performed [H] and indeed showed an- 
harmonicities. However, in the region of O-tt transition, 
the obtained CPRs are not odd functions of flux as they 
should be, indicating that the physics is spoiled by other 
effects [281 [221 • 

Here we report on the successful measurement of the 
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Figure 1: a Scanning electron microscopy image of the mea¬ 
sured asymmetric SQUID (see text and [32|). b Differential 
conductance normal state of the CNT junction 

versus gate voltage Vg and bias voltage Vsd for a large range 
of Vg. Figures c and d focus on two Kondo zones called re¬ 
spectively A and B. A IT magnetic field is applied to destroy 
superconductivity in the contacts. The number of electrons in 
the last occupied energy levels is indicated in white. In light 
blue, -^-^iysd) at zero bias is plotted (axis on the right). In 
the normal state, the reference JJ contribution is a constant 
which was subtracted to obtain the plots. The white symbols 
in c correspond to the theoretical fit of the conductance (see 
text). 


CPR of a CNT-based hybrid junction over the entire 
O-TT transition. This constitutes the first experimental 
demonstration of the O-tt transition controlled by the su¬ 
perconducting phase ip. A very important part of our 
analysis is the comparison between the measured CPRs 
and theoretical ones, computed for the Anderson model 
with superconducting leads using a numerically exact 
quantum Monte Carlo (QMC) method. The excellent 
agreement provides strong support that we indeed ob¬ 
serve the transition resulting from strong local electronic 
correlations. 

We fabricated a CNT-based QD, connected to super¬ 
conducting leads and embedded in an asymmetric mod¬ 
ified SQUID (Fig. [^a). This device, a SQUID contain¬ 
ing the QD JJ (here the CNT) and a reference JJ with 
critical current high compared to the one of the QD JJ, 
allows us to determine the CPR of interest [3I1I32]. The 
switching current Is of the SQUID versus magnetic flux 
is measured. The CPR of the QD JJ is then obtained by 
extracting the modulation of Ig around its mean value 
{Is). Our device possesses a second reference JJ and a 
third connection as described in Ref. [32]. This allows 
us to characterize each junction independently at room 
temperature, and to measure both the CPR of the CNT 
and its differential conductance in the superconducting 
state. 


The CNTs are grown by chemical vapor depo¬ 
sition on an oxidized doped silicon wafer [33] . 
A three-j unctions SQUID is constructed around 
a selected nanotube with the following materials: 
Pd(7 nm)/Nb(20 nm)/Al(40 nm), AlOx and Al(120 nm) 
mi- The sample is thermally anchored to the mixing 
chamber of a dilution refrigerator of base temperature 
50 mK and measured through low-pass filtered lines. A 
magnetic field B is applied perpendicular to the loop to 
modulate the phase difference across the CNT-junction 
by 27r5S'/4>o, with the superconducting flux quantum 
4>o = h/2e and S the loop area. 

We first characterize the sample in the normal state, 
measuring the differential conductance dl/dVsd versus 
bias voltage Vsd for various backgate voltages V^, using 
a lock-in-amplifier technique. The contacts are made of 
Pd/Nb/Al with a gap of A = 0.17 meV ± 10%, a value 
very close to the gap of A1 but considerably smaller than 
the Nb gap because of the Pd layer. A magnetic field of 
IT is needed to suppress superconductivity in these con¬ 
tacts. Even though such a magnetic field significantly af¬ 
fects the Kondo effect, the results of Fig. [^show Coulomb 
diamonds and an increase of the conductance at zero-bias 
in some diamonds, a signature of the Kondo effect. The 
four fold degeneracy, characteristics of clean carbon nan¬ 
otubes with orbital degeneracy [Ml [35] . is clearly seen 
(Fig. 0b). This allows us to determine the dot’s occu¬ 
pancy indicated on the figure. We focused on two ranges 
of gate voltages, corresponding to diamonds with odd oc¬ 
cupancy, where non-zero conductance is observed at zero 
bias, zone A (around Vg = 6.2 V) and zone B (around 
= 1.4 V) (see Figs. 1 c and d). Among all the Kondo 
ridges leading to a 7r-junction in the superconducting 
state, those show the widest extent in gate voltage of 
the O-TT transition, which makes the measurements more 
accurate. The height in Vsd of tho Coulomb diamonds 
gives the charging energy U = 3.2 meV ±10% in zone A 
and 2.5 meV ± 10% in zone B. Due to the complex inter¬ 
play of the Kondo scale ksTK and the Zeeman energy, 
that are of comparable size, F and the contact asymme¬ 
try cannot be determined directly from the experimental 
results; theoretical modeling is required. 

The junction is modeled by an Anderson impurity 
model [TIH2Q] with right {R) and left (L) BCS super¬ 
conducting leads and superconducting order parameter 
The interaction of electrons on the QD is given 
by a standard Hubbard term with charging energy U 
and the coupling of the leads to the QD is described by 
the energy independent hybridization strength We 

solve it using the numerically exact CT-INT Monte Carlo 
method m in the normal state (A = 0) in a magnetic 
field and calculate the finite temperature linear conduc¬ 
tance for different F^^/i? as a function of the dot energy 
e, defined relative to particle-hole symmetry. The ampli¬ 
tudes of the magnetic field and the charging energy are 
fixed to the experimentally determined values B = IT 
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and U = 3.2 meV of zone A m- A comparison with the 
measured conductance at zero source drain voltage (i.e. 
in equilibrium) yields a set of parameters that fit the ex¬ 
periment best. We find = 0.44meV, Tj^/Tl =4. 

We also slightly varied T to estimate the electronic tem¬ 
perature in the sample and obtain T = 150 mK. Addi¬ 
tionally, this procedure gives a reliable way of extracting 
the conversion factor a = 39 between the applied 
gate voltage and the dot on-site energy e (as done in 
|2Q|). Our best fit is displayed in Fig. [^c (white sym¬ 
bols). After reliably estimating all parameters, using the 
formula Tk = y^TU/2exp{—7r ^^^^~^ ^ ) [36] . we can di¬ 
rectly show that for zone A the dot is indeed in the regime 
of strongest competition between Kondo correlations and 
superconductivity with ksTK ~ A (Fig. [^a). 
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Figure 2: a Modulation of the switching current of the 
SQUID, proportional to the CPR of the CNT junction, ver¬ 
sus the magnetic field for various gate voltages in zone A. b 
Modulation of the switching currents near the transition for 
several gate voltages, c CPR extracted from the previous data 
and rescaled by 1.33 (see the text) near the transition (green 
continuous line). The theoretical predictions resulting from 
QMC calculations are shown as black lines. The dashed lines 
are guides to the eyes and represent the contributions of the 
singlet (0-junction, in blue) and the doublet state (7r-junction, 
in red). 


Next, superconductivity is restored by suppressing the 
IT magnetic field and the CPR is measured in both 
Kondo zones, extracting the modulation of the switch¬ 
ing current 51 s versus magnetic field (Fig. from the 
critical current of the SQUID. To measure the switching 
current, the SQUID is biased with a linearly increasing 
current with a rate ^ = 37 /iA/s and the time at which 
the SQUID switches to a dissipative state is measured. 



Figure 3: a Modulation of the switching current of the SQUID 
for zone B, which exhibits a partial transition from 0- to tt- 
junction, b CPR extracted from the previous data near the 
transition (green continuous line). The dashed lines are guides 
to the eyes and represent the contributions of the singlet (0- 
junction in blue) and the doublet state (vr-junction in red). 


This process is reproduced and averaged around 1000 
times, the whole procedure being repeated at different 
values of magnetic field below a few Gauss, small enough 
to preserve superconductivity. To obtain the modulation 
of the switching current dig versus magnetic field, the 
contribution of the reference junctions (around 90 nA) 
is subtracted. As demonstrated in Ref. [32], 51 s is pro¬ 
portional to the CPR of the CNT junction. It should 
be noted that this kind of system, in particular near the 
O-TT transition, is very sensitive to the electromagnetic 
environment, which therefore needs to be optimized [28] . 

The main results of this work are presented in Figs. 
an di where we show the extracted CPRs for gate values 
over the entire transition regime (curves c.l to c.6 of Fig. 
I^for zone A and curves b.l to b.3 of Fig. |^for zone B). 
We now analyze qualitatively the shape of these curves. 

On the edges of Kondo zone A, far from the transition 
(Fig. I C.l), the junction behaves as a regular JJ with 
a CPR proportional to sin((/ 2 ) (0-junction). In contrast, 
at the center of the Kondo zone (Fig. c.6) the CPR 
is TT-shifted {51 oc sm{(p + tt)) and has a smaller ampli¬ 
tude characteristic for a 7r-junction. In between, the CPR 
is composite with one part corresponding to 0- and an¬ 
other part to TT-junction behavior. The latter first occurs 
around cp = tt^ giving rise to a very anharmonic CPR. In 
the middle of the transition region, we find period halving 
(Fig.|c .4) (38). This evolution of the CPR between a 0- 
and TT-j unction is consistent with the finite temperature 
transition of the dot’s magnetic state (between singlet 
and doublet) which is controlled by the superconducting 
phase difference [121 Ull H] • 

A more precise analysis of the transition allows to at¬ 
tribute this O-TT transition to a competition between the 
Kondo effect and the superconductivity. Indeed, around 
the center of Kondo zone A (Fig. [^a), the 7r-junction 
extends over a range of 60 mV of gate voltage. Accord- 
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ing to the dl/dVsd in the normal state (Fig. and the 
conversion factor a, the odd diamond has a width in 
gate voltage of about 82 mV, larger than the 7r-junction 
regime. Consequently, this 0 to tt transition is not simply 
due to a change in the parity of the dot filling but to an 
increase in the ratio A/Tk (see Fig. [d] a). This is even 
more obvious for Kondo zone B (Fig.T^b) where the 0 
to TT transition is incomplete. 

For a quantitative comparison between theory and ex¬ 
periment, we performed a second CT-INT calculation in 
the superconducting state {B = 0) for zone A [37] to ob¬ 
tain the CPRs in the transition regime. We used the mea¬ 
sured value of the superconducting gap A = O.lTmeV 
and the previously determined parameters and computed 
the Josephson current as a function of the phase differ¬ 
ence ip. The theoretical CPR are calculated at various e 
(related to 4^ by e = aVg) and plotted as black lines in 
comparison to our experiments in Fig. [^c2 to c5. Since 
our setup yields a switching current that is necessarily 
smaller than the supercurrent, the experimental CPRs 
were multiplied by a unique correction factor chosen to 
obtain the best agreement with the QMC results. The 
agreement for the shape of the CPR is excellent; how¬ 
ever a shift of the energy level Se = 0.28 meV of the 
theoretical CPRs is needed to superimpose them with 
the experimental ones m- The QMC calculations pre¬ 
dict a transition region centered around a smaller e than 
measured experimentally (see supplementary materials); 
a deviation between experiment and theory which we cur¬ 
rently do not understand. Note however that the width 
of this transition is very well reproduced. 

The comparison of the measurements and calculations 
can even be refined employing Fourier decompositions 
I{(p) = ai sin((/:?) + a 2 sin(2(^) + as sin(3(/?) + ... of the 
27r-periodic CPRs. The first three amplitudes suffice to 
describe the experiment perfectly; see Fig. b where 
they are shown as functions of e. The theoretical model 
thus exactly captures the nontrivial finite temperature 
phase dependence of the measured Josephson current. 

An important information that can also be extracted 
from the experiment is the gate voltage dependence of 
the critical phase (pc at which the system switches from 
0 to TT-junction behavior, i.e. the CPR has 0-behavior 
for cp G [O^cpc] and 7r-behavior for cp G [(pc^‘^7T — (pc] 
(Fig. I c). At T = 0 this switching at cpc is asso¬ 
ciated to a first order level crossing transition and ap¬ 
pears as a jump of the current from positive to negative. 
For T > 0 the transition is washed out and the CPR 
is smoothed. However, at small enough T, cpc depends 
only weakly on T {cf. Rnfs.fTTl ITOl 130] b In Fig. c, 
we compare ipc{^) from experiment and theory. Both 
display the same characteristic shape, that can be un¬ 
derstood based on the atomic limit of the Anderson im¬ 
purity model with A ^ F. A straightforward extension 
of the T = 0 atomic limit calculation (such as presented 
in ini) to the case of asymmetric level-lead couplings 



£ (meV) 

Figure 4: a Calculated Tr (black line) is compared to A 
(dotted line) as a function of the level energy e (e = aSVg 
with 6Vg the gate value measured from half-filling and a = 
39 meV/V). In green dots, the critical phase pc is plotted, 
defined as the phase, different from 0 and tt, for which the 
CPR equals zero (see text and inset of c). The O-tt transition 
occurring when pc switches from 1 to 0, the system is indeed 
in the regime ksTR ~ A. b Fourier analysis of the CPRs 
of zone A near the transition for different level energies e. 
The circles correspond to the experimental CPRs whereas the 
continuous lines correspond to the QMC calculation (red, blue 
and black : harmonics 1, 2 and 3). c Measured pc (red circles) 
versus e. pc extracted from the QMC calculation and shifted 
is also shown (blue squares). The black line shows the result 
of a two-parameter fit of the analytical curve obtained in the 
atomic limit (see text). 

gives pc{^) = 2 arccos — (e/h)^. For the experimen¬ 
tal A ^ F the dependence of g and h on the model pa¬ 
rameters cannot be trusted; we rather fit both {g = 2.2, 
h = 0.72 meV) and obtain very good agreement. This 
shows that the e-dependence of pc is a strong character¬ 
istic of the O-TT transition. 

In conclusion we have experimentally shown that the 
magnetic state of a CNT quantum dot junction, singlet or 
doublet, can be controlled by the superconducting phase 
difference. This has to be contrasted to previously mea¬ 
sured gate controlled transition. It is achieved by probing 
with unprecedented accuracy the evolution of the Joseph¬ 
son current at the O-tt transition. We have shown that 
the CPR has a composite behavior, with a "0" and a "tt" 
component and that the phase at which this transition 
occurs is gate dependent. The measurements are success¬ 
fully compared to exact finite temperature QMC calcu¬ 
lations. The possibility to measure precisely the CPR of 
correlated systems motivates the study of systems with 
different symmetry such as Kondo SU(4) [39] or with 
strong spin-orbit coupling |4Q| . 
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SUPPLEMENTARY MATERIALS 
Sample fabrication 

The carbon nanotubes (CNT) are grown by Chemical 
Vapor Deposition (CVD) on an oxidized doped silicon 
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wafer [33]. The wafer is used as a backgate. Then the 
three-junctions SQUID (with area S ^ 40 fim^) is de¬ 
signed on top of it by electron beam lithography, with 
resin bridges for the tunnel junctions. A first trilayer of 
Pd(7 nm)/Nb(20 nm)/Al(40 nm) is deposited with an 
angle 15° (Fig. [^. Its superconducting gap is found to 
be inhomogeneous: the CNT, which is in contact with 
the Pd, sees a gap A = 0.17 meV while the gap for the 
tunnel junction is A = 0.63 meV . Then the A1 is oxi¬ 
dized and a layer of Al(120 nm) is deposited on top of it 
with an angle -15 ° . This layer has a superconducting 
gap A = 0.18 meV. 

The reference JJs obtained have tunnel resistances 
Rni = 1.6 kU and Rn 2 = 1-9 kU, from which we esti¬ 
mate the expected critical currents m- Id = 330 nA 
and Ic 2 = 280 nA. 



6 nm Palladium + 20 nm Niobium 

+ 40 nm Aluminium /-15° deposition 500nm Silicon oxyde (insulator) 

A10„ (insulator) ^ ... .... 

IWBI Doped silicon wafer (backgate) 

I I 120 nm Aluminium / 15° deposition 

Figure 5: Summary of the different layers constituting the 
sample 


Additional data for the O-tt transition 

On Fig. I^is plotted without any adjustment the super¬ 
current at a fixed superconducting phase difference (p for 
zone A, as a function of the level energy of the dot. While 
the experimental data are represented with the red line, 
the dotted red line represents QMC data. Here, only the 
positive values of e have been plotted. The O-tt transition 
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Figure 6: Comparison of raw data for the O-tt transition 


happens for a smaller e in theoretical simulation than in 
the measurement. 

Critical phase 

The transition from the 0-junction to the 7r-junction is 
a first order quantum transition and is linked to a change 
of the ground state from a doublet to a singlet state on 
the quantum dot at a critical value of the phase pc- At 
finite but low temperatures, we only see a washed out 
signal of this behavior but can nevertheless try to ex¬ 
tract the phase boundary as a function of gate voltage. 
We define the value of the critical phase at temperature 
T by the conditions I{pc,T) = 0 and ^{pc^T) < 0. 
At the particle-hole symmetric point e = 0 and at suf¬ 
ficiently low temperatures pc is known to only weakly 
dependent on temperature [T7| and will therefore provide 
an accurate estimate of the phase boundary. In order to 
check how pc depends on temperature away from the high 
symmetry point e = 0, we have performed a CT-INT cal¬ 
culation at different temperatures (expressed by inverse 
temperatures P = and for the other parameters 

as in the experiment. Our findings show that pc indeed 



Figure 7: Effect of temperature on the critical phase calcu¬ 
lated current phase relation at O-tt transition (e = O.SmeV) 
for the parameters of the experiment. The CPR has been cal¬ 
culated for different temperatures from 145 to 580 mK. For 
temperatures low enough (lower than 200 mK), the critical 
phase Pc does not depend anymore on the temperature. 

depends on temperatures at high temperature but con¬ 
verges to the zero temperature value at low ones. The 
results displayed in Fig. demonstrate clearly that at 
the temperature T = 150 mK relevant to our experiment 
(corresponding to ^ 77 ^^^) Pc essentially represents 
the ground state result. This shows that the result shown 
in the main text is in fact the accurate ground state phase 
diagram. 









































